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Abstract 

We study the fluctuation properties of the asymmetric simple ex- 
clusion process (ASEP) on an infinite one-dimensional lattice. When 
iV particles are initially situated in the negative region with a uni- 
form density p_ = 1, Johansson showed the equivalence of the current 
fluctuation of ASEP and the largest eigenvalue distribution of ran- 
dom matrices. We extend Johansson's formula and derive modified 
ensembles of random matrices, corresponding to general ASEP initial 
conditions. Taking the scaling limit, we find that a phase change of 
the asymptotic current fluctuation occurs at a critical position. 
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The one-dimensional asymmetric simple exclusion process (ASEP) is one 
of the most transparently analyzable models in non-equilibrium statistical 
mechanics PjJ 121 El H] • in the last two decades, many exact results have been 
worked out for the stationary and dynamical behavior of ASEP. Among all, 
it is now known to belong to the Kardar-Parisi- Zhang (KPZ) universality 
class 

The study of the fluctuation properties of ASEP has been recently re- 
newed after Johansson's work[6j. Johansson evaluated the current fluctua- 
tion of ASEP for an initial condition with a step at the origin, using the 
connections to random matrix theoryjTj by way of combinatorics of semis- 
tandard Young tableaux[8j. The key ingredient of Johansson's argument is 
an integral formula, which evidently shows the equivalence of the current 
fluctuation and the largest eigenvalue distribution of random matrices. 

Following Johansson's work and relying on Baik and Rains' result pi I1U| 
ITT] on combinatorics, Prahofer and Spohn|l2] made a plausible conjecture for 
a phase diagram of the fluctuation properties. In their scheme, the univer- 
sality classes of random matrices characterize the ASEP fluctuation phases. 

In this paper, we derive a generalization of Johansson's integral formula, 
which enables us to treat general initial conditions and establishes a con- 
nection to modified ensembles of random matrices. An asymptotic analysis 
of this new formula yields a phase diagram in agreement with Prahofer and 
Spohn's conjecture. 

Our approach is based on the determinantal formula for the Green func- 
tion found by Schutz[T3]. Suppose that N particles labelled 1, 2, • ■ ■ ,N from 
the left start from yi, y 2 , . . . ,yN at time 0. The probability that particles are 
at Xi, X2, ■ ■ ■ , xn at time t is given by 

P(x u x 2 , ■ ■ • , x N ; t\y u 2/2, • • • , y N ; 0) 
= det[F fc _j(x fc - yf, t)}^ k=1 , (1) 

where 

k=o v 7fe 

Here (n)& = n(n+l) ■ • ■ (n+k—1). The function F n (m, t) has some interesting 
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properties such as 

F n+1 (m + l;t) -F n+1 (m + l;0), (3) 
F n+l (m;t). (4) 

It is also remarked that F n (m,t) can be written in terms of the confluent 
hypergeometric function (or Laguerre function), though it is not explicitly 
stated in [T3] . 

The Green function is already a nontrivial result. However, in order to 
further study physical quantities, one still has to take summations over initial 
and final configurations. Our first observation is that a certain summation 
also results in a determinantal expression: 

Proposition 1. Suppose that N particles labelled 1,2, •• • ,N from the left 
start from yi, y 2 , . . . ,ijn at time 0. The probability that the particle j (1 < 
j < N) has hopped at least Xj — yj steps before time t is given by 

^2 p ( x ^ x 2, ■ ■ ■ , x N ; t\yi, 1/2, • • • , Vn; 0) 

Xj<Xj (l<j<N) 

= dtA[F h - j+1 {X k -y j ;t)]» k=1 . (5) 

This can be proved by using successively. Notice that the only difference 
between (JIJ) and (0) is the suffix of F n : in the latter it is bigger than that in 
the former by one. 

In this article we are interested in the fluctuation of the integrated current 
or the position of a particle. For this purpose, one only needs the special case 
of ((nj) with Xj — yi = M + j — 1 (1 < j < N), which we state as a corollary. 

Corollary 2. Suppose that N particles start from yi,y2, ■ ■ ■ ,Vn at time 0. 
The probability that the leftmost particle has hopped at least M steps before 
time t is given by 



,x N ;t\yx,y 2 , ••• ,yjv;0) 

M)£=i- ( 6 ) 




dsF n (m; s) = 



oo 
l=m 



Ppi > yi + M] 
= ^2 P(xi,x 2 ,--- 

yi+M<xi<X2<-<x N 

= det[F k _ j+1 ( yi - Vj + M + k - 
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Next we derive another formula for the same quantity. To obtain the 
result, we need 

Lemma 3. If /(ijv, tw-i, • • • ,ti) ^ a totally antisymmetric function of N 
variables, the following holds. 

*t N-i r t N _ j+1 N-i tN _ k 



/ (l/ v J [ / (l/ v j J J / <1-xy /,.: 

Jo j=1 Jo k=1 Jo 

k rsN-k,i-i 

X ] [ / d-SN-k,lf(tN, Sjv-1,1, SAf-2,2, • ' ' , Sl,7V-l) 

1=2 J ° 

N ~ 1 rt N ~ 1 rt N - k k rs N _ k ,i-i 

= J | / dt N -j ] i / <l>'v JJ / dsjv /,./ 

j=0 ^ fc=l ^ «=2 J ° 

x /(%) Sjv-1,1, Sat-2,2, - - - , Sl,7V-l) 
1 /•* 

=fF^-in/ d ^ n (^-w(w ••,*!). (7) 

11.7=1 .7- jr=0 ^° l<j<k<N 

Proof. Using a mathematical induction on N, we first prove a slightly 
generalized version of the first equality 

/ dt N / dt N -i ■ ■■ dh 
Jo Jo Jo 

X / dSjv-1,1 / dSAT_ 2 ,l / dsjv o2- 



/ ds A r_ lj i / dsiv_2,i / dsj 
Jo Jo Jo 

rti fsi,i 

/ d Sl;1 / 
Jo Jo 



x / dsu / dsi, 2 --- / dsi jiV -i 

X f(tN, Sjv-1,1, SjV-2,2j ' ' ' , Sl,AT-i; ^l) 

= / dt N / dtAr_i • • • / dh 
Jo Jo Jo 

' I dSiV-li / dSjV-2,1 / CISjv-2,2 



ftN-1 ftN-2 PSN-2,1 

/ dsjv-i,! / ds N _ 2 ,i / ds 
Jo Jo Jo 

rti r s i,i rsi,N-2 

x / dsi,i / dsi >2 --- / dsi jA r_ 
Jo Jo Jo 



'o Jo Jo 

X f(tN, SiV-1,1, SN-2,2, ■ ■ ■ , Sl,AT-i; tl), (8) 
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in which the function / also depends on ti. 
When N — 1, the statement is trivial. 

Now we assume that the statement is true for N — 1 and show that the 
statement is true for N. If we divide the region of the integration over t\ on 
the right hand side (RHS) as f* dti = J* 2 dt± + f* dti, the term corresponding 

to f* dti is rewritten as 

/ dt N [ d^-i ■■■ [ dt 2 [ dti 
Jo Jo Jo Jt 2 

ftN-l ftN-2 fSN-2,1 

X / dsjv-l,l / dsjv-2,l / dsAr.2,2 " " " 

Jo Jo Jo 

/ ds 2 ,i / ds 2 , 2 • • • / 
Jo Jo Jo 



r*2 f s 2,l f s 2,N-3 

x / ds 2 ,i / ds 2 , 2 • • • / ds 2 ,v-2 
'o Jo Jo 

ftl f s l,l /"Sl,AT-2 

x / dsi,i / dsi >2 --- / dsi iA r_i 
Jo Jo Jo 

X f(tN, SjV-1,1, SN-2,2, ■ ■ ■ , $2,N-2, Sl.AT-l! ^l) 

= / dt N [ dtN-!" - [ dt 3 [ dh 
Jo Jo Jo Jo 

X / dSjV-1,1 / dsjv-2,l / ds N _ 2 ,2 " " 

Jo Jo Jo 

x I dt 2 I ds 2 .i / ds 2 , 2 • • • / ds 2 ,v-2 



rti pt2 r s 2,i 

/ dt 2 / ds 2 ,i / ds 
Jo Jo Jo 

rti rsi } i rs\ t N-2 

/ dsi,i / dsi j2 • • • / dsi 
Jo Jo Jo 



x / dsi,i / dsi j2 • • • / < L.s-j \ 
'o Jo Jo 

x f(tN, Sn-1,1, Sn-2,2, ■ ■ ■ , S 2 ,N-2, S^jV-l! ^l)- (9) 



This is zero since / is antisymmetric in s 2 ,v-2 and si,at-i while they are 
integrated in the same way. Hence the region [0,t] of the integration over t\ 
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on the RHS of (jSJ) can be replaced by [0, i 2 ]- One obtains 
the RHS of © 

dt N / dtiv-i • • • / dt 2 
'o Jo Jo 

X / dSAT-1,1 / dSAr_2,l / dSAT_2,2 " " ' 

Jo Jo Jo 

x / ds 2 ,i / ds 2 , 2 • • • / ds 2 ,N~2 
Jo Jo Jo 

X #(tjV, SjV-1,1, S/V-2,2, •-• ,S2 tN - 2 ;t 2 ), (10) 

where 

g(t N , Sjv-1,1, ■ ■ • , S2,at_2; ^2) 

"*2 /"tl /" s l,l PSl,N-2 



dti / ds M / dsi >2 ■ • ■ / dsi j7V _i 
'0 ./o Jo Jo 

X /(^AT, SjV-1,1) $N-2,2, • • ■ , S2,N-2, Sl.AT-l! ^l)- (H) 

We note that <7(t/v, s;v-i,i, • • ■ , $2,n-2] h) is antisymmetric in iV — 1 vari- 
ables t N , Sjv-1,1, • • ■ , s 2 ,iv_2- From the assumption of mathematical induction, 
J Q dt 7v J dt/v-i • • - J dt 2 in (fTTH) can be replaced by J Q dtN f Q N dt^-i • ■ ■ J Q 3 dt 2 . 
Using the definition of g, we arrive at the desired statement. 

Next the second equality of the lemma is proved by a mathematical in- 
duction on N. 

When N — 1, the statement is trivial. 

Let us assume that the statement is true for N — 1 and show that the 
statement is true for N. Let us first define 

h(tN, sn-i,i, • • • , S2,jv-2J 



= dti ds M / dsi i2 • • • / dsi iA r_i 
Jo Jo Jo Jo 

X f(t N , Sjv-1,1, -5iv_2,2, ■ ■ • , Sl,jV_i). (12) 

This function is antisymmetric in iV— 1 variables t N , s^-i,i, Sjv-2,2, • • • , S2,n-2- 
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It can be further rewritten as 

h(tN, Sjv-1,1, • • • , S 2t N-2) 

dsijv-i / dti / ds u / dsi , 2 • • • / ds liiV -2 



'0 J 8l,Kt-l J S 1,N-1 J Sl,N-l J S 1,N-1 

x f(tN, S N-1,1, Sn-2,2, - - - , S2,AT-2? Sl,7V-l) 



TT7 777 / dSi,N-l(t — Sl t N-l) N 1 f {tN , S N-l,l, S N -2,2, s 2,N-2, Sl,N -l) 

{1\ — 1J! Jq 



1 



(^-l)Uo 



dti(* - ti)^ V(*JV> SiV-l,l> S AT-2,2, • • • , S 2) iV-2, (13) 



Using the assumption of mathematical induction, we rewrite the middle ex- 
pression of (J7J) as 

-t rt rt 

dt N / d%_i ■ ■ ■ / dti 
'0 Jo Jo 

'tiV— 1 ftN-2 rSN-2,1 

dSjV-1,1 / dSAr_2,l / dSAr_2, 2 ' ' ' 

'o Jo Jo 

X / ds 2i l / ds 2 ,2 • ■ ■ / ds 2 ,AT-2 

Jo Jo Jo 

x / dsi,i / dsi )2 • • ■ / dsi,jv-i 
Jo Jo Jo 

j s Ar-l,l> SiV-2,2; ■ ■ • , S2,AT_2, Sl,N-l) 
-t rt rt 

dt N / dtjv-i ■■■ dt 2 
'o Jo Jo 

"ijV-l ptN-2 rSN-2,1 

dSiV-1,1 / ds W _2,l / dSAr-2,2 • • • 

'o Jo Jo 

ft2 f s 2,l fS2,N-3 

X / ds 2i i / ds 2i2 --- / ds 2 ,7V-2 

Jo Jo Jo 

X h(t]Sf, Sn-1,1, Sn-2,2, ■ ■ ■ , S2,N-2) 

= ^jv-i , / dt N / dtjv-i • • • / dt 2 i [ (t k - tj) 

Uj=l V- ^0 J Jo 2 <j<k<N 

x h(t N ,t N -x, ■ ■ ■ ,t 2 ). (14) 
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From (JT3J), we find that it is equal to 



— ^ — - / dt* / dtjv-i • • • / d*i J] - %) (t - ti)^- 1 
ll i= iJ'!^o 7o ./o 2 < 3<fe<i v 



]— I dt N [ d^_i-.. / dt 1 ^(-l)" 1 - 1 

=1 j!^o ./o ./o m=1 



X 



n (tk-tj) {t-t m ) N - i f(t N ,t N ^--- ,t 1 ). (is) 

l<j<k<N 

Thus the second equality of the proposition is established because of the 
identity 

N 

^(-lr- 1 n (t k - 1,) (t - u^ 1 = n 

m=l l<j<k<N l<j<k<N 

j,k^m 

This identity can be proved by observing that the left hand side is antisym- 
metric polynomials in £i,i2> • • • , tjv and that the coefficient of t^ -1 on both 
sides is (-l)^ 1 n 2 <i<fc<ivOfc ~ tj)- n 
Now we state one of the main results of this paper: 



Theorem 4. 



F[xi > yi + M] 



-^n/ dt ^ n (tk-tMtN,--- M (17) 

llj=lJ- j=o ^° l<j<k<N 



where 



1 AT 



= det[F_ i+ i(j/i - + M - 1; t N „ k+1 )Yl k=1 . (18) 

When M < Un — Hi, the RHS of should be understood as its analytic 
continuation with respect to M . 
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Proof. First let us assume that the condition M > Vn ~ Vi is satisfied. 
Applying Lemma |3] to the RHS of (JBJ), we have 

^2 P ( X li x 2: • ■ ■ , %n] t\y x , y%, • • • , Vn\ 0) 

yi+M<xi<X2<---<X]si 

-t nt N rtz 

dt N / dt N -! ■■ dti 



JO JO 



x det[F fc _f(3/i - + M + k - 2; tjv-fe+i)]^ 

dt N / dtjv-i ■ ■ • / dti 
o v/o Jo 

/■*iv-i ri-N-2 rstf— 2,1 

X / dSjV-1,1 / dSAT_ 2 ,l / dsAT_ 2 ,2--- 

Jo Jo Jo 

x / dsi,i / dsi >2 ■ ■ ■ / dsi iA r_i 
Jo Jo Jo 

x f(tN, SjV-1,1, Sn-2,2, ■ ■ ■ , Si t N-i) 

i r* 

= 771^11 dtN -J II (tk-tj)f{t N ,--- (19) 

[[j=lJ- j=0 J° l<j<k<N 

Next we would like to consider the case M <y N — y x . Assume, however, 
M > yN — Vi for a moment. Let us denote by Gjsr(t) the RHS of (fT7|) with / 
replaced by fo, which is defined by 

foitN, tN-1, ■ ■ ■ ,ti) 

= det[F { yi — yj + M — 1; t N ^ k+1 )]f >k=1 . (20) 
By expanding the determinants, one obtains 



GnW = > sgn orr — ^ — — 

^ r(m -2/r(iV) +Af) 

7(yi-yr(i) + M + tr(iV)-l,t) 
T( yi - y T(1) + M) 



x 



(21) 



Here a and r are permutations of 1, 2, • • • , N. The function 7(2?, t) is the 
incomplete Gamma function defined by 



t 

s „x—l 



7(2, t) = / e-'^-Ms, (22) 
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where the real part of x is positive. Just as in the case of the ordinary 
Gamma function T(x) = j(x,oo), the incomplete Gamma function can be 
analytically continued via the contour integral representation 



7^ = ^ e-^dz, (23) 
e 1 JC(t) 



where C(t) is the contour starting from t, enclosing the origin anticlockwise 
and ending at t. This expression is meaningful for all x on the complex 
plane except at the locations of poles. In fact it is known that 7(2, t) has an 
expression, 

7(M) = E 1 x V (24) 
nUx + n) 

n=0 v ' 

and hence has precisely the same singularities as T(x); the poles are located 
at x = 0, —1, —2, • • • and the corresponding residues are 1, —1, 1/2, • • • . Tak- 
ing these facts into consideration, one realizes that the RHS of (}2*Tj) can be 
analytically continued to the region M < y N — yi, because each (incomplete 
or ordinary) gamma function can be analytically continued and the singu- 
larities of the gamma functions on the numerators are cancelled by those on 
the denominators. 

By using (jljl. one can rewrite / as a linear combination of fo and vice 
versa. Hence the RHS of (fT7|) can also be analytically continued to the 
region M < yN — Hi- The value for an integer M is finite and gives the 
desired probability F[xi > y\ + M]. □ 

For the special case y d = -N + j (j = 1, 2, • ■ ■ , N), the RHS of (fTTj) 
reduces to 



N-l „t N 

rM-N-t 



J M,N f = L Q Jo .J^KN 7=1 



e-*0, (25) 



where Zm,n = Ylf=ijK^ ~ N + j ; — 1)!. This is the Johansson's integral 
formulajn], the RHS of which yields the largest eigenvalue distribution of the 
Laguerre unitary ensemble of random matrices. One can therefore analyze 
the ASEP current fluctuation by means of random matrix theory. Johans- 
son derived this formula by way of combinatorics of the semistandard Young 
tableaux. The derivation of our theorem is physically more transparent, be- 
cause we have just summed up the Green function. Furthermore our theorem 
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generalizes Johansson's formula and allows us to study arbitrary initial con- 
ditions. In addition, the RHS of (|17|) can be considered as a modification of 
matrix integrals appearing in the random matrix theory [Jj. 

The initial condition in the above theorem has a fixed configuration of 
N particles. One can also consider a superposition of initial conditions by 
summing up over particle configurations. In fact, it can be even technically 
easier to work with a superposition. Notice that Gjv(t) defined in the proof 
of the above theorem represents the probability that the leftmost particle has 
hopped at least M steps before time t for a superposition of initial conditions. 
The technical advantage of replacing / by fo is that it can be written in a 
form 



When M < yjy — yi, the RHS should be again understood as special cases of 
its analytic continuation with respect to M. This point will be discussed in 
some detail afterwards. 

To proceed further, we construct monic polynomials Cj(t) and Dj(t), so 
that Cj{t) (Dj(t)) is a linear combination of Ck{t) (dk(t)) with k — 0, 1, • • • , j 
and satisfy 



(C 3 , D k ) = / C 3 (t)D k (t)w(t)dt = h 3 5 jk , j, k = 0, 1, 2, ■ • • , N - 1, (27) 



Jo 

where hj (j — 0, 1, • • • , N — 1) is a constant. After some reflection, one sees 
that Cj and Dj can be expressed as 




N-l 



N 



(26) 



Cj(t) oc 



co(t) 

*(t) 



(c , d ) 
(ci,d ) 



(c , dj-x) 



(28) 



(cj,d ) 




Dj{t) oc 



d {t) 
di{t) 



(do, c ) 
(cfi,c ) 



(do, Cj_i) 
(di, 



(29) 



dj{t) (dj,c ) ■■■ (d^Cj-x) 

with the condition that they are monic (the coefficients of Cj(t) and dj(t) are 
!)• 
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Let us now define 

( j )] (t) = ^=C j {t)^Mt), (30) 



^(t) = ^=D ] {t)^Mf), (31) 
'Al- 



so that (|27j) reads 

/■oo 

4>j{t)i) k {t)dt = 5 jk . (32) 

Then we can show that Crjv(i) has an expression 

G N (t) = det(l- K N ), (33) 

where the determinant on the RHS is the Fredholm determinant and the 
kernel K N (x,y)xi(y) is given by 

N-l 

K N (x,y) = J2 ( l ) j(^(y), (34) 

3=0 

X*W = { 0, otherwise (35) 

and I = (t, oo). This Fredholm determinant expression can be derived by 
following the strategy of Gaudin[l4j. 

In the following, we consider a special choice 

yi = - N+ l-Y, y J = -N + j (j=2,3,---,iV) (36) 

in (|2Up . where Y is a nonnegative integer. This choice is one of the simplest 
generalizations of the case studied by Johansson. The corresponding initial 
condition is that the particles are distributed in the negative region with a 
uniform density 

and that the particle density in the positive region is zero. As we noted, the 
fluctuation was already conjectured by Prahofer and SpohnfTJ]. Our formula 
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Cj(t) 


= t j , 


djit) 


I* ' 


wit) 


= t a e- 1 



allows us to prove their conjecture. In terms of Cj(t) and djit), the choice 
(J2HJ) corresponds to 

= 0,1,---,AT-1), (38) 

,- = o,i,...,jr-2), 

= iV-l, 

(40) 

with i(T = iV + y — 1. For the case specified by (jSHD"©, C 5 {t) is the monic 
version of the Laguerre polynomial, 

Cjit) = kJ x Lf\t), ij = 0, 1, • • • , N - 1), (41) 

with kj = (—iy/j\. In terms of Cj(t)'s, the orthogonality of the Laguerre 
polynomials reads 

POO 

/ Cjit)C k it)wit)dt = hf ) S jik , (42) 
Jo 

with hf ] = j\T(M-K+j). The polynomial Dj(t) can also be equated with 
the monic Laguerre polynomial for j = 0, ■ ■ ■ , N — 2, 

Djit) = kJ l Lf\t), ij = 0, 1, • • • , N - 2). (43) 
Hence, for j = 0, 1, ■ • • , N — 2, is identical to /i^: 

^ = hf = j\F{M -K + j), ij = 0, • • • , N - 2). (44) 

As for Djv_i(t), one uses the determinant formula (|2*9j) . After some compu- 
tation, one obtains 

!>„-.(«) = wt ^J^/ ft («) 

The explicit formulas for Cjv-i(t) and D^-iit) allow one to compute fr/v-i 
in (|2*7)l . The result is 

h N -i = (-l) JV - 1 r(M)(-K)^_ 1 . (46) 
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As we mentioned below ()26p . for the case M < y N — y 1 , some care should 
be taken for the interpretation of the formulas appearing in the above dis- 
cussion: the orthogonality relation (|27jl has its meaning as its analytic con- 
tinuation. For j — 0, 1, • • • ,N — 1 and k — 0, 1, • • • ,N — 2, (J2T)) is equivalent 
to the orthogonality relation of the Laguerre polynomials, 

fJ^+a+Tj £ j a) (') i **'(*)*"e"'dt = (47) 



Using a formula 



rM m _f (-lfflj + q + l) g 

^r(j-/ + i)r( a + / + i)r(/ + i)' 1 8J 



one finds 

r(j + 1; 



r(j + a + 1 



/■oo 

/ Lf\t)L^\t)t^dt 
Jo 



(-l) m+T T(j + l)r(fc + a + l)r(m + n + a + 1) 

As a function of a, there is no singularity on the RHS. 

For j — 0, 1, • • • , N — 1 and = iV — 1, on the other hand, one has to 
use the expression f!45|) of -Djv-i. However, since only the j = iV — 1 term 
in ()45|) gives a nonzero contribution, the orthogonality relation (J2T)) can be 
analytically continued in a similar manner. 

Finally we consider the scaling limit of the current fluctuation for the 
special case (j36|) with a uniform initial density p_ in the negative region. 
The results are given in terms of the transcendental functions J~2(s) and 
.7-1 (s), which were originally introduced in random matrix theory. The largest 
eigenvalue distribution of Gaussian (or Laguerre) unitary and orthogonal 
ensembles are written in terms of J~2{ s ) and J-\{s), respectively. 

If one introduces the function u(x) which satisfies the Painleve II equa- 
tion, 

d 2 

-u = 2u 3 + xu, (50) 



with the asymptotics 



dx 2 

u(x) ~ Ai(x) x — > oo, (51) 
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the functions ^(s) and J-\(s) are represented as [TBI ITo] 



^2(5) = exp 
Ti(s) = exp 



(x — s)u(x) dx 



2 



x — s)u{x) 2 dx 



u(x)dx 



(52) 
(53) 



They also have the Fredholm determinant representations fT7\ ITSj 

F 2 {s) = det(l-AT 2 ), 



(J- 1 ( S )) 2 = det(l-K 1 ). 
The kernels K2(^,rj)xi(v) an d v)Xi(v) j respectively, are given by 

Ai(flAi'fa) - Aifa)Ai'(0 



(54) 
(55) 



/■oo 

iir 1 (e^) = -K'Airy(e^) + Ai(0 / Ai(77-A)dA 

./0 



(56) 
(57) 



and J = (s, 00). 

Now let us consider the asymptotics of the current fluctuation. We set 
M = Y + 7^V (7 > 0). This corresponds to measuring the current at a 
position a = (7 — 1)N. The results depend on the initial density on the 
negative sites, p_. Now the results of our asymptotic analysis (in agreement 
with Prahofer and Spohn's conjecture^]) are summarized as 

Theorem 5. Set t = cN + dN^-^s where 

c=(l + v ^) 2 , d = T l, \l + VlY l?J ■ (58) 
(i)For UmN^ooY/N < ^y, Gpf{t) tends to F 2 {s) in the limit N — > 00; 

lim G N (t) =F 2 {s). (59) 

N—*0O 

(u)When Wm^Y/N = ^ (i.e. p_ = 1/(1 + a = N{1 - 2p_)/p 2 _), 
Gpf(t) tends to Ti{s) 2 in the limit N — > 00/ 



lim G N (t) = (^i(s)) 2 . 



(60) 
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Proof. Let us rewrite K N (x,y) as 



JV-l 



K N {x,y) = <t>f\x)<t>j (v) + 4>n-i(x)iI>n-i(v) 



where (f>f\x) = (f>j(x)(j = 0, • ■ ■ , N-2), ^(x) = C N ^(x) y/Mxj/^h^. 
Using the Christoffel-Darboux formula for Cj, the first term is rewritten as 



JjM Y)N ^IZ^l^l VN wwn-W) (62) 

x-y 

Here and in the following we need the asymptotic behavior of the Laguerre 
polynomials. In terms of the function defined by 

= (-1)" fr4m V2 V^)L { n a) (x), (63) 
\{a + n)\J 

one has 

^v.W-^iAiK-A) (64) 

(x = cN + dN 1 ' 3 ^), which can be obtained by applying the saddle point 
method to the integral representation of <pffl{x). Then it is not difficult to 
see that the first term in ()61jl multiplied by dN 1 ^ 3 tends to K^ iTy (^,ri). 

We next consider the second term on the RHS of (|61|). One sees that, 
when liniTv^oo Y/N < y/j, the term with j = N — 1 becomes dominant in the 
limit N — > oo. Using this expression and the asymptotics of (p in (|64p. we 
can show that the second and the third terms in (|61|) are negligible. Hence 
the case (i) is proved. 

When Y = one uses the asymptotics, (JM|), of ip to find 

poo 

dJV^Vjv-iO^-iG/) ~ Ai(0 / dAAi(r? - A). (65) 

The third term in (JfiHj) is again negligible. This leads us to K±(£, rf) in (|5T|). 
ending the proof of the case (ii). □ 
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The density profile of ASEP in the scaling limit t — > oo (with N = t/c) 
is known to bej^j 



We can therefore conclude that the current fluctuation is described by ^(s) 
in the region with a linearly decreasing density. At the singular position x c = 
(1 — 2p-)t, a phase change occurs (a Gaussian fluctuation is conjectured [T2] 
at x < (1 - 2p_)t). 

Let us finally remark that, owing to the particle- hole symmetry of ASEP, 
we can extend the result to a measurement at —a. In that case, the initial 
particle density in the negative region is fixed to 1. 

To summarize, we have studied the fluctuation properties of the asymmet- 
ric simple exclusion process (ASEP) on an infinite one-dimensional lattice. 
We first obtained a determinantal formula by summing up the Green func- 
tion. This can be turned into a form of a modified random matrix ensemble. 
As an application of our formula, we considered the fluctuation of the in- 
tegrated current when the particles are initially distributed in the negative 
region with a uniform density p_. A generalization of the technique in ran- 
dom matrix theory has been utilized to study the asymptotic properties of 
the model. Consequently we found a phase change of the current fluctuation 
at a critical position. 

The authors would like to thank Prof. Makoto Katori, Prof. Hideki 
Tanemura and Mr. Takashi Imamura for useful discussions and comments. 
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